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Abstract. - We present results from our simulations of biopolymer translocation in a solvent which explain 
the main experimental findings. The forced translocation can be described by simple force balance arguments 
for the relevant range of pore potentials in experiments and biological systems. Scaling of translocation time 
with polymer length varies with pore force and friction. Hydrodynamics affects this scaling and significantly 
reduces translocation times. 
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The transport of biopolymers through a nano-scale pore in 
a membrane is a ubiquitous process in biology. For example, 
in protein import into mitochondria, chloroplasts, and perox- 
isomes the translocation occurs with the aid of a membrane 
potential [ 1 ] . Experimental work on forced (or biased) translo- 
cation is largely motivated by finding methods for reading the 
DNA and RNA sequences. These nanopores are typically ei- 
ther fabricated solid-state [2, 3] or a-hemolysin (a— HL) pores 
in lipid bi-layer membranes [4, 5], Foundation for the the- 
oretical work was laid in the classic treatment by Sung and 
Park [6], which was based on the assumption that the poly- 
mer segments on the two sides of the membrane reside close 
to separate thermal equilibria. However, the validity of this 
approach was questioned already in [7, 8], where the authors 
noted that the pore force regime in which the polymer's relax- 
ation time towards equilibrium is smaller than the characteristic 
translocation time is marginal and that the approach would be 
invalid even in the unforced translocation for sufficiently long 
polymers. Theoretical work, inconsistent with experiments, has 
since evolved in different directions. 

The role of computer simulations has been largely to support 
the theoretical work which neglects hydrodynamics. Hence, 
results from simulations where hydrodynamic interactions are 
included are few and, due to their being computationally de- 
manding, often fairly qualitative [9, 10]. In addition, the gener- 
ally used Monte Carlo method gives unphysical behaviour for 
larger pore force values relevant for experiments and biological 
systems [2,4,5], as we have shown [11]. Very recently multi- 
scale simulations on biopolymer translocation in a solvent were 
reported to give results in accordance with experiments [12, 13]. 



Our motivation for the present study is two-fold. First, by 
using realistic dynamics we want to find explanation for the dy- 
namics of the experimentally observed translocation processes. 
Secondly, we want to determine the effect of hydrodynamics 
on forced polymer translocation, previously studied only in the 
unforced case [9, 10]. We use a hybrid multi-scale method, 
where the polymer follows detailed molecular dynamics and 
the coarse-grained solvent stochastic rotation dynamics (SRD). 
The solvent is divided into cells, within which fictitious solvent 
particles perform simplified dynamics where collisions among 
them and with the polymer beads are taken effectively into ac- 
count by performing random rotations of the random part of 
their velocities, v l {t + At SRD ) = H[vi(t) - v cm (t)] + v cm (t), 
where v.- L are the particle velocities inside a cell, AAsrd is 
the time step for solvent dynamics, R is the rotation ma- 
trix, and v cm is the centre-of-mass velocity of the particles 
within the cell. Hydrodynamic modes are supported over the 
cells. Optionally, they can be switched off by not adding back 
v cm after the random rotation, which is particularly feasible 
for pinning down the effect of hydrodynamics. The above- 
described collision step is followed by the free-streaming step 
n(t + M SRD ) = n(t) + Vi(t)At SRD . Thermostating is done 
by rescaling all solvent particle velocities so that equipartition 
theorem is fulfilled at all times. More detailed descriptions of 
the method can be found e.g. in [14-16]. 

In this paper we study the forced translocation where the two 
sides separated by walls are not hydrodynamically coupled. To 
achieve this we use a non-aqueous pore, i.e. there are no sol- 
vent particles inside the pore. This corresponds closely to the 
experiments we aim to model and also addresses the theoreti- 
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cal predictions, where the two subspaces separated by the wall 
are taken to be uncoupled. In addition, technically speaking 
the coarse-grained solvent dynamics does not allow for overly 
confined spaces, but implementing solvent dynamics in dimen- 
sions smaller than the SRD cell dimension would require de- 
tailed molecular dynamics for the solvent. The linear SRD cell 
dimension in our model is Ax = 1.0 = b, where b is the poly- 
mer bond length. 

The standard bead-spring chain is used as a coarse-grained 
polymer model [11, 17]. Adjacent monomers are connected 
with anharmonic springs, described by the finitely extensible 
nonlinear elastic (FENE) potential, 
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Here r is the length of an effective bond, R = 1.5 the maximum 
bond length. The Lennard-Jones (LJ) potential 
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0, r > 2- 1/6 a, 



(2) 



is used between all beads. The parameter values were chosen 
to be e = 1.2, a = 1.0, and K = 60/er 2 . The used LJ potential 
with no attractive part mimics good solvent condition for the 
polymer. 

We do not include a harmonic bending potential, which 
would change the above-described freely-jointed chain (FJC) 
to the worm-like-chain (WLC). Elastically, the FJC and WLC 
were seen to differ only marginally and their hydrodynamic 
characteristics were found identical in the present model [17]. 
The swelling exponents measured for the FJC and WLC in the 
present model were identical, which means that since packag- 
ing of polymers is not an issue in the present geometry, as e.g. 
in capsids [18], the bending potential does not constitute an im- 
portant factor to the translocation process. 

The model geometry contains a slit 1 formed by two walls 
perpendicular to the x-direction. Periodic boundary conditions 
are applied in y- and z-directions. A third wall, three polymer 
segments, b, thick and impermeable to the solvent, is placed in 
the middle of the system, perpendicular to the z direction. No 
slip boundary conditions are applied between the three walls 
and the solvent. A circular nanopore of diameter 1.2 & is placed 
in the centre of the middle wall. The force / acting on the 
beads inside the nanopore is constant and local for the pore, 
which models well the experimental setups and biological sys- 
tems, where solvents are good ionic conductors eliminating any 
potential gradients outside the pore. The polymer beads inside 
the pore are not coupled with hydrodynamic modes or with the 
heat bath and in the directions perpendicular to the cylindri- 
cal pore walls experience a damped harmonic potential Uh, de- 



scribed by -VUh = F h 



v , where k = 1000, 



c = 10, t x y is the polymer bead position with respect to the 
centre axis of the cylindrical pore, and v XiV is the velocity com- 
ponent perpendicular to the pore walls. Thus Uh centres the 



'The slit mimics experimental setups including confining walls. The slit is 
large enough not to affect dynamics of the translocating polymer. 



polymer along the z-directional axis of the pore. The potential 
is chosen large, so no hairpin configuration can enter the pore 
as its width is effectively small. Hence, the polymer segment 
inside the pore remains rather straight. In the z direction, the 
polymer beads experience either zero or finite friction in the 
pore. In the zero friction case the polymer beads inside the 
pore are moved by a constant force included in the molecular 
dynamics and the momentum is conserved in the SRD step per- 
formed every 50th MD step. In order to investigate the effect 
of pore friction the momentum was optionally destroyed in the 
SRD step in a manner similar to switching off hydrodynamic 
modes in the solvent. This amounts to the polymer beads expe- 
riencing finite friction inside the pore. 

Diffusion of the chains and single particles has been verified 
to satisfy equipartition theorem in equilibrium. The swelling 
exponent value of a self-avoiding-chain, v = 0.6 ± 0.05, was 
measured and the fluctuations of the radius of gyration the poly- 
mer configurations in directions aligned with and perpendicular 
to the side walls were measured and found equal [17]. Thus the 
side walls do not affect the dynamics of polymers. Translo- 
cation simulations were started from initial configurations that 
were checked to be in equilibrium with respect to the radius 
of gyration, R g . The measured i? 9 's were clearly smaller than 
the channel dimensions, which excludes effects due to polymer 
confinement. Also, finite (simulation box) size effects were 
verified not to affect the dynamics by measuring the relaxation 
time T r ~ (R g (t)R g (0)) for polymers of length N = 200 
and various volumes (box sizes). The relaxation times are of 
the same order than the largest translocation times presented in 
Fig. 1 b). The simulation box sizes are [25, 32, 32] for N < 50, 
[32, 32, 32] for N < 100, [40, 32, 32] for N < 200, [40, 40, 40] 
for N < 400, and [60, 60, 60] for N = 800. The measured 
radii of gyration are 4.91 ± 0.01 for N = 50, 7.19 ± 0.01 for 
N = 100, and 10.93 ± 0.01 for N = 200. 

In order to characterise the native translocation process, we 
use pore force values that are sufficient to induce translocation, 
hence obviating any additional constraints for preventing the 
polymer from sliding back to the cis side (see a snapshot of 
a translocating polymer in Fig. 2 c)). As has been noted, ad- 
ditional constraints can potentially change the observed scal- 
ing of the translocation time r with polymer length N [19]. 
We compare our simulated results to experiments [2, 4, 5] and 
also to those obtained analytically from the Brownian translo- 
cation framework [6, 20] and the numerical results supporting 
anomalous diffusion [8, 10,21]. The assumption of the whole 
or parts (called folds) of the polymer close to the pore being in 
equilibrium is crucial to the application of the aforementioned 
frameworks. We show that these assumptions are invalid in 
the experimental pore force range and that hydrodynamics has 
a profound effect on the forced translocation. This, in turn, 
we show to be a highly non-equilibrium process governed by 
a simple force balance closely related to the one presented by 
Storm et ah, albeit with important modifications. 

We present the results from our simulations in reduced, di- 
mensionless units. The unit of length is defined as the polymer 
bond length, b, which corresponds roughly to the Kuhn length 
of the translocating polymer in Si-units, b. For our freely- 
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Fig. 1: (Colour online) a) The distribution of translocation times r for chains of length N = 100 and a constant pore force / = 10, with 
hydrodynamics (A) (averaged over 300 runs) and without hydrodynamics (■) (averaged over 400 runs). Inset: Average translocation time r 
as a function of the driving force /. The scaling r ~ f a is obtained with (■) a = -0.940 ± 0.013 for / > 3, and (A) a = -0.994 ± 0.008 
f° r / > 1 without and with hydrodynamics, respectively. The chain length is N = 100, and the pore is frictional. b) & c) Average 
translocation time r as a function of the chain length TV with b) frictional pore and c) frictionless pore. Results are displayed both with and 
without hydrodynamics (HD). b) Shown are r ~ N 13 for forces / = 1, 10, 100 (from top to bottom). For forces / = 1,3, 10, 40, 100, 
(3= 1.25±0.02, 1.26±0.02, 1.31±0.02, 1.43±0.03, 1.52±0.03 without HD and f3 = 1.05 ±0.02, 1.04 ±0.01, 1.07±0.01, 1.13 ±0.01, 
1.18 ± 0.02 with HD. c) Shown are r ~ N 13 for forces / = 1, 10 (from top to bottom). For forces / = 1, 3, 10, 100, (3 = 1.50 ± 0.04, 
1.50 ± 0.03, 1.58 ± 0.03, 1.70 ± 0.03 without HD, and (3 = 1.23 ± 0.03, 1.26 ± 0.02, 1.33 ± 0.02, 1.48 ± 0.02 with HD. 



jointed chain (FJC) the Kuhn length can be taken as b = 2X p , 
where A p is the persistence length, roug hly 40 A for a single- 
stranded (ss) and 500 A for a double-stranded (ds) DNA [22]. 
In the simulations the force is exerted on three beads resid- 
ing simultaneously inside the pore. The pore force per bead 
in Si-units, /, is obtained from the dimensionless force per 
bead, /, as / = fksT/b. The simulations were performed 
at ksT = 1, which we take to correspond to T = 300 K. 
Hence, the dimensionless force / = 1 corresponds to the total 
pore force f to t = 3/, giving f tot ps 0.12 pN for ds DNA and 
hot ~ 1.6 pN for ssDNA. A typical experimentally used po- 
tential driving a polymer through the pore for both the ssDNA 
in the a— HL and dsDNA in the solid state pore is ~ 120 mV, 
which would give a pore force of ~ 50 pN for ssDNA and 
~ 110 pN for dsDNA. When charge reduction due to Manning 
condensation is taken into account, the effective force for ds- 
DNA in the solid-state pore was evaluated to be in the range 
20 - 50 pN [2, 23]. For ssDNA in an a-HL pore the charge 
reduction was evaluated to be even more drastic due to con- 
finement in the pore in addition to the normal charge reduction, 
giving / ~ 5 pN [4, 24]. This would suggest that also the esti- 
mated force for dsDNA could be smaller. The translocation of a 
polymer across a pore in a biological membrane involves in ad- 
dition friction and e.g. interaction of the polymer with the pore 
proteins [25,26], which without detailed information on those 
interactions makes exact mapping of the pore force values used 
in the simulation to those in experiments impossible. The pri- 
mary control parameter is the total pore force, f to t = Mf, 
where M is the number of points on the polymer contour on 
which the pore force, /, is exerted. On dsDNA these points 
can be taken to reside at intervals determined by the nucleotide 
spacing, which is 3.4 A for dsDNA and ps 4 A for ssDNA. 
The pore force per bead in the experiments may be estimated 
as / = zq*V/L, where the pore potential V = 120mV and 
the number of elementary charges e per nucleotide is z — 2 for 
dsDNA and z — 1 for ssDNA. The effective charge q* is taken 



as e for dsDNA [2] and O.le for ssDNA due to charge reduc- 
tion [24]. This gives / ps 1.92pN for dsDNA and / ps 0.37pN 
for ssDNA. Since the length of the solid state pore is 20 nm, 
M ps 59, ftot ~ H3 pN for dsDNA, but could be considerably 
smaller as pointed out above [2]. The length of the a-HL pore 
is 52 A, so M ps 13, giving f fot ps 5 pN for ssDNA [4]. So in 
summary, in spite of the intricacies involved in estimating the 
true force exerted on the polymer inside the pore, the experi- 
mental force magnitudes are included in the pore force range 
/ £ [1, 100] used in our simulations and, what is more impor- 
tant with respect to the observed of out-of-equilibrium effects 
pertinent to the forced translocation, the minimum pore force 
used in our simulations is well under the minimum pore force 
magnitudes used in experiments. 

For the frictional pore we obtained the translocation proba- 
bility of P tr = 0.12 ± 0.05 for / = 0.25 that can be taken 
as a crude estimate for the minimum force required for generic 
forced translocation. P tr reached unity at / ps 1, For a fric- 
tionless pore P tr ~ 0.5 with / ~ 1, and P tr — 1 at / ps 2. Ex- 
periments on protein translocation across inner mitochondrial 
membrane, where the pore is highly frictional, showed satu- 
ration of Ptr with the pore potential 80 mV [25]. The poten- 
tial in related experiments varied from 150 to 240 mV [27, 28]. 
Hence, in terms of generic translocation, the pore force used in 
our simulations is in a range relevant for experiments. The sim- 
ulated pores were non-aqueous. Having solvent inside the pore 
may change the minimum pore force value inducing translo- 
cation. This will be investigated in a future publication. Ex- 
perimentally, an average velocity of (v) ~ 1 cm/s was mea- 
sured for a dsDNA translocating across a solid-state pore with 
a typical potential of ~ 120 mV (corresponding to a pore force 
of 20 — 50 pN) [2, 3]. We obtained an average velocity of 
v = 0.004 for / = 1, so the dimensionless simulation time unit 
corresponds roughly to i ~ 0.1 /xs. In our translocation simula- 
tions SRD step is performed every 50th MD step, so the solvent 
and polymer time steps are AtsRD = 0.1 At and AtMD = 
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0.002Ai, respectively. We obtain r\ k, 15.776 for the viscosity 
of our model [29]. For a polymer of length TV = 100 we mea- 
sured radius of gyration of R g fa 7. 198. Hence we calculate for 
the Zimm relaxation time [30], i.e. the time it takes for the en- 
tire polymer to relax to an entropically and sterically favourable 
configuration, t z = 0.398 nR 3 g /k B T w 2.4 • 10 3 . Estimat- 
ing the friction coefficient for our model [29] to be ( « 25, 
we obtain for the corresponding (Rouse) relaxation time with- 
out hydrodynamics [30] t r = ((Nb) 2 /3ir 2 k B T « 8.5 • 10 3 . 
So, in our model the estimated ratio of relaxation times with 
and without hydrodynamics of a polymer of length TV = 100 
is R — t z /t r w 3.5. We measured R w 1.6. The single- 
particle Reynold's number for this velocity in our model [29] is 
Re = 0.005, which is in the relevant regime for physiological 
solvents. 

First we determine the translocation time, r, as a function 
of pore force, /. For the scaling exponents a defined as 
t ~ f a we obtain a = -0.940 ± 0.013 for / e [3, 100] and 
a = -0.994 ± 0.008 for / e [1, 100] without and with hy- 
drodynamics, respectively, see the inset of Fig. 1 a). Hence, 
essentially r ~ / _1 was obtained, as was to be expected for 
force values large compared with thermal fluctuations. 

The distribution of translocation times, tt(t), for polymers 
of length TV = 100 is shown in the main part of Fig. 1 a). Due 
to the larger polymer velocities in forced translocation the ef- 
fect of hydrodynamics on forced translocation is much more 
pronounced than what has been seen with unforced transloca- 
tion [9, 10]. We obtain a reduction in translocation times due 
to hydrodynamics, which was also seen by Fyta et al [13, 31]. 
In addition, hydrodynamics not only significantly speeds up 
forced translocation but also reduces the variance of measured 
translocation times, which is induced by the long range corre- 
lations due to hydrodynamics, mediating the effect of the pore 
force along the polymer contour. 

The measured translocation times scale with polymer length, 
r oc TV^, both with and without hydrodynamics, see Figs. 1 b) 
and c). It is noteworthy, however, that there exists no single 
scaling, but (3 varies with pore force, /. In our simulations, (3 
starts from unity and increases with /. At constant /, smaller 
(3 was obtained for the frictional pore. To distinguish between 
the change of (3 due to increasing translocation velocity, v, and 
due to frictional term, when / was increased, scaling of t with 
TV for a pore with no friction was measured and it was found 
that the change of (3 was still significant, see Fig. 1. The exper- 
imentally obtained (3 w 1.27 for a solid state pore [2] would be 
obtained in our model with a pore force / ~ 3. Hence, it can 
be concluded that the change of (3 with / arises not only from 
the change in the frictional contribution to the translocation dy- 
namics, but also because of dynamic changes due to the change 
in v, which is a clear indication of out-of-equilibrium effects. 
In comparison, Fyta et al. [13] obtained (3 = 1.28 ± 0.01, and 
(3 = 1.36±0.03 with and without hydrodynamics, respectively, 
for the pore force / = 1, which closely corresponds to / = 1 
in our simulations as the pore length in [13] is approximately 
36. A pore of very low friction was used in these lattice Boltz- 
mann (LB) simulations. Accordingly, the obtained scaling ex- 
ponent is in fair agreement with ones we have obtained for the 



frictionless pore with hydrodynamics. Also the increase of (3 
when hydrodynamics is switched off qualitatively agrees with 
our results. Fyta et al. report the scaling exponent only for a 
single pore force magnitude. The verification of possibly non- 
universal scaling exponents varying with the pore force using 
the lattice Boltzmann method would be most valuable. 

Using Langevin dynamics (where hydrodynamics is ex- 
cluded) we checked that linear scaling (J3 = 1) can be achieved 
for any constant pore force with large enough friction. Lin- 
ear scaling has been seen with a-HL pores, whose diameter is 
smaller and friction larger than those of solid state pores. The 
important notion is that (3 varies with pore friction even at mod- 
erate friction values, regardless of hydrodynamic interactions. 

Most of our results are for a frictional pore, which is the 
more realistic. We use pore force magnitudes / > 1 at which 
translocation takes place with very high probability (P tr ~ 1) 
and hence we do not address unforced translocation. With- 
out hydrodynamics, the translocation with a frictionless pore 
approaches the scaling r ~ N 1+u w TV 16 as / increases. 
f3 = 1 + v is the scaling exponent predicted by the Brownian 
translocation framework, independent of /. (3 can be increased 
further by applying an unrealistically large pore force. This in- 
crease is due to crowding on the trans side, discussed further 
below. 

When hydrodynamics is allowed, the polymer segments are 
moved from their initial equilibrium positions already before 
actually being pulled by the tightening polymer contour. This 
is seen in Fig. 2 a), where the squared distance, R 2 e (n), of 
the polymer bead, labelled n, measured from the pore on the 
cis side as a function of the number of translocated beads, s, is 
shown. In the absence of hydrodynamics, the segments towards 
the free end are seen to remain immobile until they are pulled 
towards the pore, whereas due to hydrodynamic interactions the 
distance of the labelled bead n from the pore is seen to start de- 
creasing right from the beginning of the translocation. Hence, 
the initial configuration shows less in the translocation, when 
hydrodynamics is included. Instead, the configuration on the 
cis side continually changes towards increasingly extended one. 
Regarding only the dynamics on the cis side, for an initially 
completely extended polymer asymptotically (3 — > 1 as / is in- 
creased, which explains the reduction of (3 at constant / when 
hydrodynamics is applied. Evidently, mechanisms like squeez- 
ing of equilibrated folds through the pore, suggested in [21,32], 
contradict with this mode of motion, where alignment of seg- 
ments in the vicinity of the pore is to be expected and, indeed, 
seen in the snapshots, see Fig. 2 c). 

In order to identify the underlying mechanism in the forced 
translocation, we extract from the measured distances of the la- 
belled polymer beads from the pore, R pe the number of mobile 
beads, s m . We define a labelled bead as mobile if this mea- 
sured distance, averaged over several runs, changes apprecia- 
bly. The number of mobile beads, s m , as a function of translo- 
cated beads s is read off from the inflection points in Fig. 2 a) 
depicting the measured R 2 e (s/N) in the case of no hydrody- 
namic interactions. In Fig. 2 b) s m is plotted as a function 
of translocated beads s when hydrodynamics is not included. 
Linear dependence s m = ks is obtained. Up to lengths of 
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Fig. 2: (Colour online) a) Averaged squared distances of beads num- 
bered 30, 50, 70, and 90 from the pore as a function of the number 
of translocated beads s for polymers of length N = 100 with and 
without hydrodynamics, b) The number of mobile beads, s m , (see 
text) vs number of translocated beads, s, both normalised to the poly- 
mer length, TV, in the case of no hydrodynamics. / = 10 in a) and 
b). c) 3D snapshot of a translocating polymer of length N = 100 at 
s = 35. Here the pore is frictionless and / = 2. For clarity reasons, 
walls are not shown, d) The radius of gyration R g on the trans side 
as a function of the number of translocated segments s. Chain length 
TV = 100. The applied pore forces for the curves from top to bottom 
are 1, 3, 10, 20, 40, and 100. Measured equilibrium R g 's for s = 50 
and s = 100 (■). The pore is frictional in a), b) and d). 



N ph 200, k ~ N~ x and levels off to a constant value > 1 for 
longer polymers. At all times, the drag force, fd, balances with 
the constant pore force, fd is exerted on mobile beads, so in the 
absence of hydrodynamics fd ~ s m (v), where (v) is the av- 
erage velocity of the mobile beads. When the whole chain has 
translocated, fd ~ N m (v), where N m — kN. With no hydro- 
dynamics, the beads are set in motion from their equilibrium 
positions, so the distance d of the mobile bead furthest from 
the pore scales as d ~ N v . The average translocation time then 
scales as r - (d)/{v) ~ kN 1+v ~ N 1+u ~*. For the data 
in Fig. 2 b), where / = 3, we obtain x ~ 0.35 that accords 
with the measured (3 = 1.26, see Fig. 1 b). With the pore force 
/ = 100 the fc's for s m = ks are smaller and the measured s m - 
s curves for different N appear more aligned. Asymptotically, 
k — ► 1, ViV, as / — > oo, i.e. polymer beads are translocated at 
the same rate that they are set in motion. Removing the friction 
from the pore also makes k values smaller and more identical 
for different TV due to translocation becoming faster. Both the 
increase in the pore force and reduction in the pore friction take 
the scaling exponent (3 toward 1 + v due to s m and hence the 
drag force, fd, remaining more constant throughout the translo- 
cation. Hydrodynamics changes the form of the drag force, fd 
no more depends strictly linearly on s m for configurations of 
moving polymer segments, but all beads are set in motion in 
the beginning of translocation. Yet, the above described mech- 
anism is still clearly underlying the translocation also when hy- 



drodynamics is included, see Fig. 2 a). For a constant / hydro- 
dynamics reduces (3 due to enhancing collective motion of the 
polymer towards the pore. Hence, the initial equilibrium posi- 
tions do not determine j3 like in the absence of hydrodynamics. 

In addition to the above described dynamics on the cis side, 
there is a potential contribution from the crowding of the poly- 
mer beads close to the pore on the trans side. This indeed 
can be seen from the snapshot in Fig. 2 c) and the measured 
radii of gyration, R g , of the translocated parts of the polymer 
on the trans side, which are clearly smaller than the corre- 
sponding equilibrium R g even for the smallest force values, see 
Fig. 2 d). The crowding is more enhanced for longer polymers 
thus increasing [3, in agreement with our measured transloca- 
tion times, Figs. 1 b) and c). Asymptotically, for N oo, 
the average velocity of the polymer beads on the cis side, (v), 
would have to diminish as s rn increases to maintain the force 
balance fd=fp = constant. Eventually s m would be so large 
that the polymer barely moves, which would bring the translo- 
cation to the regime where diffusive motion of the beads shows. 
For finite polymers this is not the case. Instead, the simulated 
polymers whose lengths measured in Kuhn lengths are well in 
the range of polymers used in experiments clearly show that 
the moving polymer segment does not slow down to velocities 
where diffusive motion could be seen. Crowding was also re- 
ported in the LB simulations in [13]. However, the alignment 
of segments in the vicinity of the pore, Fig. 2 c), was not ob- 
served, but on the contrary polymers were reported to stick to 
the wall on the cis side. 

In conclusion, we have studied forced polymer translocation 
by a model where hydrodynamics is taken judiciously into ac- 
count. In our minimal model the pore is non-aqueous, which 
precludes hydrodynamic coupling of the two chambers sepa- 
rated by the wall, and the pore potential is the only driving force 
for the translocation. No additional mechanisms for preventing 
the translocated polymer segments from sliding back to the cis 
side were included. Using this model, the smallest pore force at 
which polymers translocate was estimated to be in the order of 
1 pN, according with experimental findings [2,3,25]. The used 
pore force values cover the biologically relevant range and thus 
characterise well the essential dynamics of forced translocation 
in biological systems, and DNA experiments. Hydrodynamics 
was shown to significantly speed up the translocation and di- 
minish variation in the translocation times. The scaling expo- 
nent, (3, of the translocation time with respect to polymer length 
was seen to increase with pore force, /. The obtained scaling 
exponents and their variation with the pore force could be ex- 
plained by simple force balance at the pore and the observation 
that the rate, at which the size of the part of the polymer in 
cis side moving toward the pore grows with respect to the part 
translocated to the trans side, varies with /. A simple estimate 
was given for the case when hydrodynamics is not included, 
which agreed with the numerical results. Consequently, a sin- 
gle universal exponent cannot describe translocation for all /. 
The magnitude of pore friction was also seen to change (3. Lin- 
ear scaling, (3 = 1, was obtained for large enough friction. (3 
was shown to change also due to crowding mechanism on the 
trans side. As crowding increases with pore force it is an ad- 
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ditional mechanism accounting for the increase of (3 with /. 
In summary, by using realistic dynamics where hydrodynam- 
ics is included we have shown that experimentally observed 
forced translocation can be described by a simple force bal- 
ance. The forced translocation process was also shown to be a 
highly non-equilibrium process for the experimentally relevant 
force regime, which explains the discrepancy between theoret- 
ical approaches and experiments. We have shown that no uni- 
versal scaling of the translocation time with the polymer length 
exists. Experimental verification of this by using different pore 
potential magnitudes would be very important. 
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